Spin injection, accumulation and spin precession in a mesoscopic non-magnetic metal 

island 
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We experimentally study spin accumulation in an aluminium island with all dimensions smaller 
than the spin relaxation length, so that the spin imbalance throughout the island is uniform. Elec- 
trical injection and detection of the spin accumulation are carried out in a four terminal geometry 
by means of four cobalt electrodes connected to the island through tunnel barriers. We model the 
system theoretically and we investigate the role of the ferromagnetic electrodes on the spin accumu- 
lation at the limit at which the electron diffusion time can be neglected. We present measurements 
of spin accumulation at room temperature and at 4.2 K: in both cases the spin accumulation signal 
is larger than the Ohmic resistance of the aluminium island. From magnetisation precession mea- 
surements at room temperature, we extract a spin relaxation time Ts/ = 60 ps and a polarisation 
P = 8% for tunnel barriers with resistances as low as 20 . We show that the precession 

measurements are invariant under the interchange of voltage and current electrodes, and under the 
reversal of magnetic fields and magnetisations, according to the reciprocity theorem. We show that 
spin accumulation and spin precession in a system with uniform magnetisation can be described in 
terms of the (relative) orientation of the ferromagnetic contacts' magnetisations and we determine 
from precession measurements the angles between the magnetisation direction of the contacts. 
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PACS numbers: 85.75.-d, 73.23.-b 



I. INTRODUCTION 



Creating and manipulating a non-equilibrium magneti- 
sation in a non-magnetic metal is a central requirement 
in the field of spintronicsi. The orientation of an electron 
spin injected in a non-magnetic metal is the result of the 
interaction of the spin intrinsic magnetic moment with 
the magnetic fields in which the spin is moving. In the 
presence of a uniform magnetic field, the spin precesses 
coherently around the field's direction and its orientation 
changes with a uniform precession frequency. 

In a diffusive metal, a (non-uniform) effective magnetic 
field arises from the relativistic motion of the spin in the 
electric field of the metal ions and the defects^. This is 
called spin orbit interaction and it is responsible for the 
randomisation of the spin orientation: the relaxation of 
the non-equilibrium magnetisation occurs by transferring 
the spin angular momentum to the metal lattice in a time 
scale in the order of 100 ps. 

To induce a spin current and a spin accumulation, a 
spin dependent scattering is required for the conduction 
electrons. The usual approach in all-electrical transport 
experiments is to drive a current from a ferromagnet 
(FM) whose band structure is spin dependent, to a non- 
magnetic (N) metal. It has been noted that the main 
obstacle for efficient spin injection in the diffusive regime 
is the short spin relaxation length in the ferromagnet, a 
problem known by the name of conductance mismatch"^. 
In GMR experiments'^, with vertical devices, the useful 
signal can be made large enough for practical applications 
by reducing the distance between the ferromagnetic lay- 
ers. In lateral structures, this is not a feasible solution, 
due to technological limitations. Also, in the clean con- 
tact regime, measures have to be taken to ensure that 



either the current path is perpendicular to the FM/N in- 
terface or that the voltage probes are separate from the 
current path. Anomalous magneto-resistance and Hall 
effect can mimic and hide the spin accumulation^. 

In order to overcome these problems, tunnel barriers 
were proposed^ at the interface between the ferromag- 
netic metal and the normal metal layer (the tunnel bar- 
rier conductance being proportional to the FM density of 
states) , thereby making the tunnel barrier the dominant 
(spin dependent) resistance of the system. 

Previous experiments have studied the spin current in 
systems larger than the spin relaxation length. A semi- 
nal experiment was performed by Johnson and Silsbee® in 
the clean contact regime and four-terminal configuration, 
in a device with two lateral dimensions larger than the 
spin relaxation length. The experiment was performed 
on single crystal aluminium bar. The long spin relax- 
ation length they found, A^/ = 50 /im at 4.2 K, allowed 
them to observe a weak (tens of pV) spin precession sig- 
nal at macroscopic scale. In diffusive metallic systems, 
with typical relaxation lengths in the /im range, observa- 
tions were done by Jedema et alm/^ at room temperature 
in a one-dimensional device, both with clean contacts 
and with tunnel barriers at the FM/N interfaces in four- 
terminal devices. The spin signal in the clean case was 
about 1 mfi and about two orders of magnitude larger for 
the devices with tunnel barriers, proving their efficiency 
as spin injector/detector. Spin accumulation occurs in 
two terminal pillar structures with all dimensions shorter 
than the spin relaxation length, used to study the mag- 
netisation reversal of a thin FM layer, driven by a spin 
polarised current created by a second massive FM layeA 
The torque exerted on the FM layer is proportional to the 
spin accumulation, which is in the mV range. 
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Recently* we have performed electrical injection and 
detection of spin accumulation in an aluminium island 
with all lateral dimensions shorter than the spin relax- 
ation length As/ = yjDTgf {~ 0.6 /im at RT), where D is 
the diffusion constant and Tsf the spin relaxation time. 
Tunnel barriers separated the island from the cobalt elec- 
trodes. As opposed to the previous four-terminal exper- 
iments, the spins in our case are confined to the island, 
and since the diffusion time Tdiff — / D [L being the 
island's size) is shorter than r^/, the induced magneti- 
sation behaves uniformly within the island, so that the 
spatial variation of the magnetisation can be disregarded 
and the system is zero-dimensional with respect to the 
spin. 

The description of coherent spin transport in a diffusive 
metal is obtained from the Boltzmann transport equa- 
tion. The two-channel (spin up and spin down) model of 
Valet and FerliiS, which successfully describes the exper- 
imental results of giant magneto-resistance, is however 
limited to the coUinear (magnetisation either parallel or 
anti-parallel) case. In the general case, one has to retain 
all the information about the spin direction inside the 
bulk metal and at the FM/N interfaces. 

A theoretical approach to systematically studying the 
transport through FM/N interfaces in the non-coUinear 
situation was developed by Brataas et al}^ . The relevant 
parameter, alongside the interface conductance for spin 
up and spin down electrons , G-'- , is a (complex valued) 
mixing conductance term G^-^, describing the reflection 
of electron spins perpendicular to the magnetisation of 
the ferromagnet. G'^'^ is related to the amount of angular 
moment that the electron spin has transferred to the fer- 
romagnet, and plays an important role in the description 
of the spin torque and spin pumping^^'^'^. 

Here we present a systematic study of non-collinear 
spin accumulation in a small metallic island, extending 
the results of Ref. Theoretically, we apply the cir- 
cuit theory to the system and show how the presence of 
the FM contacts provides an additional (and anisotropic) 
mechanism for the relaxation of the spin accumulation, 
with the relaxation occurring at a faster rate in the di- 
rection perpendicular to the FM magnetisation axis. 

In Sec. we derive a formula in the zero dimensional 
limit {Tdiff ^ Tsf) for the spin accumulation as a func- 
tion of the contacts properties. Sec. IIIII describes the 
sample fabrication and Sec. IIVI is a short summary of 
the relevant theory. The geometry of the sample and the 
measuring configurations are described in Sec. [V] 

Sec. I VII presents an extensive set of measurements of 
spin accumulation in a non-magnetic island at 4.2 K and 
at room temperature (RT). We use spin precession as a 
tool to analyse the spin accumulation and extract the rel- 
evant parameters such as the spin relaxation time, and 
the direction of the magnetisation of the FM contacts. 
We also show that, for our device, the magnetisation re- 
laxation is independent of the mixing conductance term. 



II. THEORY 

To completely characterise the electronic transport in 
the linear regime, including the spin, four chemical po- 
tentials are necessary, one spin independent (charge) and 
three spin dependent (magnetisation). The main idea 
of the finite element theory of Brataas et alm^, is to 
divide the system into (normal or ferromagnetic) nodes 
connected to each other or to reservoirs by contacts (in- 
terfaces). A contact can be spin selective, that is it can 
have different conductances for the two spin populations. 

In each node, spin accumulation appears as a result 
of the spin currents through the contacts. In turn, the 
amount of spin current through each contact is deter- 
mined by its conductance and by the chemical potentials 
of the two nodes on each side of the contact. 

Thus, the problem of the transport in the system is 
broken down into the solution of the motion of charge 
and magnetisation inside a node with the additional 
boundary conditions given by the charge and spin cur- 
rent through the contacts. The finite element theory of 
Brataas et al. provides an elegant way to describe the 
charge and spin currents through FM/N interfaces. 

In the following, we will briefly review the elements 
of the theory that are relevant for our experimental sit- 
uation. We will then present an analytical solution for 
the situation in which the electron diffusion time inside 
the island Tdiff can be neglected when compared to the 
spin relaxation time Tsf. This last assumption is equiv- 
alent to saying that the distribution of magnetisation is 
uniform in the island and we will therefore call it, a zero- 
dimensional system from here on. 

In a FM/N device, and for arbitrary configuration of 
the magnetic reservoirs, a spin current injected from 
a FM reservoir to a node of normal metal will cause 
a non equilibrium accumulation of magnetisation. In 
the linear regime, we describe the transport properties 
by means of four chemical potentials ii{x),^{x), where 
fi{x) = //(e)de and H — J f(e)de are the spin indepen- 
dent and spin dependent chemical potentials in the node 
and /(e) and f(e) are the spin independent and spin de- 
pendent distribution functions, fi = {nx^ Hy, jJLz) repre- 
sents the spin accumulation in different directions and 
its magnitude. In equilibrium, no spin accumulation 
exists in normal and ferromagnetic nodes or reservoirs. 

Only in the particular case in which the spin accu- 
mulation has the same direction throughout the entire 
system (and the contacts are all coUinear), one can use a 
description in terms of two, spin up and spin down, chem- 
ical potentials: these are related to /i, /x by ^| = -I- 
and ^1 = /i - 

In a reservoir, the spin independent chemical potential 
is set by the applied bias voltage eV . 

In the experimental device, see Fig. ^ the node is 
an aluminium island 400 nmx400 nmx30 nm, the four 
cobalt electrodes act as reservoirs and the AI2O3 tunnel 
barriers at the aluminium/cobalt interface are the con- 
tacts. The motion of charge and magnetisation in the 



3 




anti-parallel). Gj and Gj are the conductances for the 

up and down spin channels, G, = G] + G'^ is the total 

contacts' conductance and Pj = (G] - Gj)/(G] + G^) 
the polarisation of the interface. 

Also the assumption that spin transport through the 
contacts can be specified only in terms of r"™, r"™ im- 
plies that a spin up electron have zero probability of being 
converted into a spin down electron, that is no spin flips 
occur inside the contacts^ ^. 

It is also assumed that spin accumulation in the FM 
side can only be coUinear to the magnetisation direction, 
i.e. the spin dependent chemical potential is of the form 
fi^ = |/j^|m, as the large exchange field rapidly ran- 
domises the spin component perpendicular to m. 

The charge current entering the normal metal readsi^ 

e2^, = G, {^lf P,G, (m, ■ /. - |/.f |) (1) 

and the spin current 



e^ij =PjGjipf ~ M)mj + Gj {fif - (m^ • /x)mj) 

+2ReGymj x {mj x /x) - 2/TOG]^mj x /x (2) 

fij being the spin independent chemical potential of 
the FM electrode j. The conductances are defined ac- 
cording to the Landauer-Biittiker formalism: 



FIG. 1: (top) A schematic representation of the system: the 
node (the square island of non-magnetic metal) is charac- 
terised by the spin independent and dependent chemical po- 
tentials /i, fi, connected to ferromagnetic electrodes, irij is 
the unit vector parallel to the magnetisation of the FM. The 
shadowed regions represent the contacts, separating the node 
from the reservoirs, (bottom) Scanning electron microscope 
micrograph of the real device. The square aluminium island is 
connected to four cobalt electrodes through transparent tun- 
nel barriers—. 



island is diffusive as the mean free path is of the order of 
5 nm at RT and 20 nm at 4.2 K. 

The finite element theory specifies the charge and spin 
particle currents i and i through the contacts that con- 
nect the island. This is related to the amplitude proba- 
bilities r"™ for the reflection of a spin up electron from 
mode m to mode n in the normal metal, evaluated at the 
normal side of the contact. If, as in our case, the tunnel 
barrier is non-magnetic and rrij is the unit vector rep- 
resenting the magnetisation direction of the j electrode, 
the natural choice for the quantisation axis is coUinear 
to nij (and f means, for instance, parallel to nij and | 



(3) 



for the spin up conductance and for the mixing con- 
ductance g'^-^ 



gt^ 



M - 



(4) 



where M the total number of modes. The spin mix- 
ing conductance affects only the component of the spin 
accumulation perpendicular to the the electrode's mag- 
netisation by rotating the spins around it, see the last 
two terms of eq. The validity of the above expres- 
sions is restricted to the case in which the contacts limit 
the total conductance^. 

In our experiment, we use FM electrodes both for the 
injection of a spin polarised current and for the detection 
of the spin accumulation. We now derive from the cir- 
cuit theory some relationships relevant in the two cases: 
FM electrode as voltage probe and FM electrode as spin 
source. 

When a FM electrode is used as a voltage probe, the 
spin independent chemical potential on the FM side is 
raised above the N chemical potential by an amount that 
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depends on the spin accumulation on the two sides of the 
contact. To see this, we set i = in eq. ^ and obtain 



where the term on the right hand side is the source 
term v and T is a 3 x 3 matrix operator 



Pj being the "efficiency" of the detector. 



(5) 



When the FM electrode is used as a spin injector, the 
charge current carries along a spin current. We use eq. 
in eq. (O, as we control in our experiment the charge 
current, and we find and expression relating ij{x) and 
ij{x): 



-2ImGYm.j x fj. + 2ReGYmj x {mj x /x). (6) 

The first term shows that a charge current I — ei car- 
ries a spin current I = /i^i, /is being the Bohr's magne- 
tone, with an efficiency given by the polarisation P of the 
interface. The second term describes a decrease in con- 
ductance because one spin channel is partially blocked. 

If the contacts have much higher resistance than the 
ferromagnetic region in which spin accumulation occurs, 
the particle currents are thus determined by the large 
voltage drop across the interface and the small spin accu- 
mulation in the ferromagnet can be neglected altogether, 
fi^ — 0. This is valid in the limit G ^ (Tf^f^ (in 
the order of Ohms for thin FM layers), where ctf and 
Af ~ 50 nm are respectively the conductivity and the 
spin relaxation length of the ferromagnet. 

For the spin independent chemical potential, one has 
to solve the diffusion equation 

-ZJVV = (7) 

with the boundary condition set by the charge current 
(its direction is along the gradient of the chemical poten- 
tial) 

ijix) = iyDD\Vfj,\ (8) 

For the spin dependent chemical potential, in the limit 
Tdiff ^ Tsf, one can neglect the diffusion term —DV^fj, 
and assume a uniform spin accumulation fia throughout 
the island. In the steady state, the injection of magneti- 
sation has to compensate for the relaxation: 
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where the first term on the right hand side describes 
spin relaxation and the second spin precession in a uni- 
form external magnetic fieldii. Using the expression for 
the spin current eq. (01, we rearrange the terms, to show 
that the presence of the contacts introduces an extra 
mechanism for spin relaxation. 

The spin dependent chemical potential can be written 
in the following formi^: 



vdV , 



(10) 



Gj{l-P^ ' 
VDe^V 



mj X fla 



^2ReGY -Gjil- Pf) 

- 2^ -> X (nij X /X,) (11) 



An explanation of the above now follows. The first 
term proportional to fia relaxes the magnetisation via 
two different mechanisms: (a) the interaction of the spin 
with the normal metal (spin orbit scattering), occur- 
ring at a rate t~j^ , and (b) the leaking of the spins to 
the leads, proportional to the interfaces' conductance, 
Gj. The time associated with the latter is the spin 

escape time Tesc = J^j^ji^ ~ P])I^Dev ■ The total 
spin relaxation time is the sum of the two contributions: 

'rel ~ ' sf ^ 'esc- 

This timescale Tesc is relevant in a two-terminal GMR- 
type of measurement: a spin dependent resistance ap- 
pears if electrons cross the second FM/N interface while 
still retaining the information about the magnetisation 
of the first FM/N interface. This is equivalent to having 

The second term plays the role of an effective magnetic 
field (jj = gfisB/h + (jJmix'- the magnetisation perpen- 
dicular to ijj precesses with constant Larmor frequency 
|a;|/27r around it. The presence of the leads introduces 
an extra term uimix, that depends on the orientation of 
the contacts and changes sign if all the magnetisations 
m.j are reversed. 

The last term affects only the spin accumulation per- 
pendicular to m^-. If this is the case, it simply becomes 
proportional to /.ta and it adds up to the spin relaxation. 
With other words, it is responsible for the anisotropic 
relaxation of the magnetisation. 
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that the coefficients 



[meGY -^Gjil-PfYj /vne^V are larger than 0, 
resulting in an enhancement of the relaxation of the 
component perpendicular to the electrodes' magnetisa- 
tion. This follows from the assumption that in the FM, 
the spin accumulation is only coUinear to the electrode's 
magnetisation. has the units of time and represents 
the spin mixing time Tmix- 

Eq. H10|) can be solved by inverting the matrix T, giving 
for ^la 
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a'^v + {lj ■ v)iA} - acj X V + J2^ CiCm, • uj)mi x v - a^,- c^nii x (nii x v) + i Y.i,j CiCj[{mi x rcij) ■ v](mj x nij) 

a3 + a|a;|2 - cja^ + (m; • a;)^] + f X), j |m, x m^P - g j fc c,CjCfe|mi • (m^- x mfe)|2 

, '_ 02) 

I 



In the experiments reported here, we use precession 
measurements as a tool to study the spin accumulation, 
by applying a uniform magnetic field to the island, for in- 
stance in the z direction. We now plot the dependence of 
the detected spin related contribution eAV = fi^' — ^J, = 

fi ■ d with d — Pjirij — Pjimj> (scaled by {lyjyeV)'^) as 
a function of the external magnetic field B = |B|z, for 
a four-contact device and the magnetic configuration de- 
picted in the inset of Fig.[2Ia) (with the magnetisation of 
the largest electrodes pointing inwards by 20° and lying 
completely on the substrate). We use renormalised pa- 
rameters g ~ G/ve^V, g^^ — G^^ jve^V and we assume 
P = 50% and Tsf = 60 ps. 

To show the dependence of eq. H12|) on g^^-, we plot the 
detected signal for arbitrarily chosen values of the pa- 
rameters g^^ and g. Figure|2Ia) shows the magnetic con- 
tribution to the total signal using g'^'^ = (1.5 -f i0.1)g, for 
different values of the interface conductance g. The width 
of the curves increases from the top trace (g = 0) to the 
bottom one {g — 16), reflecting the fact that with increas- 
ing interface conductance, more relaxation takes place in 
the leads. Fig. |2Ib) shows the signal for fixed g = 1 and 
different mixing conductancesSfi: 2(7^^ = (1.5-|-i0.1)A:, with 
k taking the values 1.5, 5 and 20. Here again the traces 
broaden, but now the relaxation of the spin is truly due 
to the mixing term 2g^^ , the third term in eq. (|ll|l . The 
solid lines represent the precession field in the z direction 
and dashed lines for the field along x. Wc also note that 
the maximum of the solid curves shifts to negative fields 
as a result of the intrinsic precession field proportional to 
2g^''-. This is not the case for the dashed traces, because 
the direction of the contacts' magnetisation leads to to- 
tal cancellation of the intrinsic precession field of the x 
component. 

The detected precession signal, also in the presence of 
the mixing conductance terms, can always be expressed 
as the sum of an "absorptive" (even) term 



and a "dispersive" (odd) onai^ 



III. SAMPLE FABRICATION 



The system under study consists of an aluminium is- 
land with lateral dimensions of 400 nm x400 nmx30 nm. 
Four cobalt electrodes of different width are connected 
to the island through tunnel barriers. A typical device 
is shown in Fig. ^ Devices are fabricated by suspended 
shadow mask technique and by electron beam (e-beam) 
lithography. We begin with a trilaycr consisting of 1.6 /im 
copolymer PMMA/MMA, 40 nm germanium and 200 nm 
PMMA deposited in this order on a 500 nm thermally ox- 
idised Si substrate. After e-beam exposure and develop- 
ment of the trilayer— , the germanium mask is suspended 
1.6 /im above the substrate. Then 30 nm Al are deposited 
under an angle to form the island, in an e-beam evapora- 
tion system with base pressure of 10^^ mbar. We notice 
that changing the evaporation rate from 0.1 nm/sec to 
0.3 nm/sec reduces the Al resistivity by a factor of 2. For 
aluminium deposited at 0.3 nm/sec, paooK / P4K — 2, and 
for 0.2 nm/sec, Pzqqk/pak = 1-3. In the following, all 
the devices have been deposited at a rate of 0.3 nm/sec, 
unless indicated otherwise. Next, we oxidise the Al in 
0.02-0.2 mbar pure oxygen for 2-5 min to produce tun- 
nel barriers (20-500 il/im^). The devices are produced 
with different tunnel barrier transparencies (from 1 kf2 
to 40 kO). Co leads (40 nm thick) are subsequently de- 
posited under a different angle. 

The devices are fabricated with decreasing tunnel bar- 
rier resistance to determine the lowest transparency for 
which the tunnel barriers still retain a sizeable spin selec- 
tivity. We also started off with the idea of measuring the 
mixing conductance term. In order to measure 
or |a;,nja;|^^ have to be comparable to the relaxation time 
Trel- It was shown" that 2ReG'T^ - G > 0, the equality 
holding true for tunnel barriers. Also, the imaginary part 
ImG^^ for tunnel barriers, is of the same order of mag- 
nitude as G. The devices we fabricated with the highest 
transparencies, have tunnel barriers of G^^ = 1 kfl and 
show a spin relaxation time of Tg/ = 60 ps and an escape 
time of Tesc = 10"^ Ts/: such a system is unsuitable for a 
measurement of G^^. 



with ^ b + SlOz and a suitable choice of Slo^ and 
. Figure |2Ic) shows the absorptive and dispersive 
terms, plotted for Scuz = 0. 



Tunnel barriers with resistances 3 orders of magnitude 
lower could not be fabricated (and probably can not) in 
aluminium oxide. The alternative would be to decrease 
the island volume by a factor 1000, but this is not feasible 
with this fabrication technology. 
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FIG. 2: Calculated spin signal contribution, proportional to 
/X ■ d as a function of the applied magnetic field |B| for a mag- 
netic configuration represented in the inset of the top panel, 
with the largest electrodes tilted inwards by 20°. In (a), gj 
are assumed to be the same for all junctions, and equal 0, 1, 
2, 4, 8 and 16 in units of voe^V (from top curve to bottom). 
For all curves P — 0.5 and Tsf — 60 ps. For the mixing con- 
ductance, we have arbitrarily chosen 2g''^ — (l.S + iO.l)^. (b) 
shows the calculated signal for a fixed value oi g = 1, now 
varying the mixing conductance 2g^^ = (1.5 + i0.1)fc, for k = 
1.5, 5 and 20. Solid lines correspond to the field applied in 
the z direction, dashed lines for the field applied along x. (c) 
The "absorptive" and "dispersive" components— as defined 
by eq. (1131 and eq. 11411 : in precession measurements, the de- 
tected spin signal can be written as a linear combination of 
these two terms, with a possible shift in the magnetic field. 



other two electrodes. For the devices with the most 
transparent tunnel barriers we could fabricate, the mix- 
ing term accounts for a correction to the spin accumu- 
lation of lO^'^. For this reason, we set Gj = gY = 
in eq. (|12|l . The following equations are derived from it, 
after some algebraic manipulation. The spin dependent 
contribution to the total signal eVs — fi ■ d: 



I 



(15) 



in the absence of magnetic field a; = 0, where s = 
Pimi — P2m2 is the source term and d = Pama — P^m/i 
the detector, if current is sent from Col to Co2 and 
voltage is detected between Co3 and Co4. vd = 2.4 x 
10^* eV^^m^'^ is the aluminium density of states at the 
Fermi energy and V the volume of the island. 

In the presence of a magnetic field, the magnetisa- 
tion perpendicular to the field precesses. When using 
ferromagnetic electrodes, one has also to be concerned 
with the stray fields generated by the electrodes them- 
selves, that thread the island. To account for these, we 
add an extra term uigt to the external magnetic field 
g^B|B|z/fi = z, a; = 6z + tj^t- 

In general 



R. = 



II -^11 



s_L • - (s^ X d±_)-u)Tsf 



1 



\'^?r^sf 



(16) 



where 1 1 means in the same direction as the exter- 
nal field. This is the z direction if also uigt has only 
a z component, uj^t.z- As mentioned before, eq. Hlt)|) 
can be written as a linear combination of the absorp- 
tive and dispersive terms (and a constant term), eq. 113|) 
and eq. 114|l . shifted in the precession field by u}st,zi with 



T_7 = r^/ + w?. 



-2 
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Eq. Hlt)|) is not invariant under the reversal of the elec- 
trodes' magnetisation s ^ — s and d — > — d, because the 
reversal also produces a change of sign uSst —^st 7^ 0. 
It obeys, however, the reciprocity relational that requires 
the interchange of voltage and current probes and the 
reversal of all magnetic fields and magnetisations. In 
the following we show that we obtain experimentally the 
same spin signal but only in the constraints given by the 
reciprocity theorem. 



MEASUREMENT CONFIGURATIONS 



IV. SPIN INJECTION EXPERIMENTS 

We first review and simplify the expressions that are 
relevant for our devices. 

Spin accumulation measurements are done in a four- 
terminal geometry: we drive a current / into and out of 
two electrodes and we detect the voltage V using the 



The four cobalt electrodes have different widths, one 
pair 100 nm wide, the other 500 nm, with the latter 
having the lowest coercive field. The magnetic shape 
anisotropy holds the electrodes' magnetisation in the sub- 
strate plane, and by applying an in-plane external mag- 
netic field along the electrodes' direction (the y direc- 
tion), we can independently reverse the direction of the 
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FIG. 3: (a) The three magnetic configurations and (b) the 
three possible independent measuring configurations: Current 
is sent from 7"*" to /~, the detected voltage is V = — V~. 
(c) The chemical potentials /i| = fi + \fi\ and /j.^ — fi — 
inside the island for the side configuration in the antiparallel, 
parallel and anomalous cases, assuming coUinear magnetisa- 
tion. We remind that only for antiparallel injectors, is 
uniform (see main text). The lines represent the spin up and 
down chemical potentials and the thick line the average 
fj,. The black dots indicate the potential measured by the 
and V~ probes for P = 1. 



magnetisation of the electrodes. We identify an antipar- 
allel configuration, in which two electrodes are pointing 
in the same direction and two in the opposite, and a par- 
allel, in which all four have the same direction, as shown 
in Fig.|2Ia). Here parallel and antiparallel are used as a 
practical shorthand notation: we will show in fact that 
the electrodes' magnetisations are non-coUinear. In the 
anomalous configuration, three electrodes (the two wide 
and one of the narrow ones) are pointing in the same di- 
rection and the fourth narrow electrode is in the opposite 
direction. 

Figure OJb) shows the three independent electrical 
measuring configurations. The current / is sent between 
/+ and /~ and the detected voltage is F = — V~ . 
The plotted signal is R = V/ 1. In the side configuration, 
the background signal is the island's Ohmic resistance. 
In the diagonal configuration, little Ohmic contribution 
is expected, owing to the symmetric position of the volt- 
age contacts with respect to the current path. The spin 
dependent contribution in the two cases is however equal 



if the island is zero-dimensional. On the other hand, the 
opposite configuration should show small spin signal as 
the widest electrodes switch at the same time (s w 0) 
and so do the narrow ones (d « 0). 

We write the total signal as the sum of a spin in- 
dependent (Ohmic) and a spin dependent contribution, 
R = Rohm+Rs- Rohm IS the island four-terminal Ohmic 
resistance and we assume it to be independent of the 
magnetic arrangement of the electrodes (we exclude for 
instance the Hall effect). 

Rs is the spin-related part. We refer to Fig. |3Jc) to 
illustrate its contribution in the three different magnetic 
configurations. Suppose for the moment that all elec- 
trodes are collinear and the barrier polarisations Pj are 
all equal. and /~ are the current electrodes, V'^ and 
V~ the voltage probes and the black dots are the volt- 
ages that would be detected if the polarisation was P = 1. 
The position of the dot on the ^| or /i| lines depends on 
the orientation of the detector. In the antiparallel con- 
figuration, a uniform non-equilibrium magnetisation in 
the island is created and — fJ.i)/2 = ^ 0. This 
potential difference is detected at the voltage electrodes 
and the signal Rs is given by eq. lfTB|) . 

In the parallel configuration, there is no net spin ac- 
cumulation. However, a spin current |I| = PI^b/s. is 
injected at and extracted at /~, giving rise to a space 
dependent magnetisation, |I| = —{(JnIJ'b/s.) ■ V(/i| — 
/i|)/2, (Jn being the Ohmic conductance of the island 
in much the same way a charge current generates a space 
dependent chemical potential, / = — (crjv/e)V/Lt. Recall- 
ing that the device is in the parallel configuration, the 
detected spin related contribution is P times the differ- 
ence of the spin up chemical potential at the and 
positions. The spin signal is a fraction of the Ohmic re- 
sistance, Rs = P^Rohm- For P = 1, the total resistance 
doubles because only one spin channel is used, halving 
the island's conductance. 

In the anomalous configuration, the V'^ probe mea- 
sures /i| and V~ detects /j,^. Owing to the symmetric 
position of V'^ and V~ with respect to /+ and /~, both 
probes measure the same amount of magnetisation and 
the spin dependent contribution is cancelled^^, i?^ = 0. 
In the anomalous configuration, one therefore expects to 
have the lowest signal, equal to the island's Ohmic resis- 
tance. 

Standard lock-in techniques are employed, with exci- 
tation currents ranging from 5 /xA to 100 /xA and with 
modulation frequencies between 4-10 Hz. 



VI. EXPERIMENTAL RESULTS 

A spin valve experiment is a four terminal resistance 
measurement in one of the three possible configurations, 
see Fig.l^fb), as a function of the in-plane magnetic field 
(in the y direction). Spin valve measurements were per- 
formed both at 4.2 K and at RT, and precession mea- 
surements only at RT. We measured 9 devices in total, 
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FIG. 4: Normalised differential resistance of a 1.8 kQ tunnel 
barrier at 4.2 K as function of the dc voltage: a peak appears 
at zero-bias and the curve shows an asymmetry. Positive volt- 
age means current flowing from Co to Al. 



one at 4.2 K only, two both at RT and 4.2 K, and 6 only 
at RT. They all show consistent behaviour. 

We report here a complete set of measurements on a 
device with tunnel barriers of i?i = 1.5 kfl, R2 — 0.90 kfi, 
i?3 = 1.6 kn and R4 = 0.75 kO. at RT^^. From these, we 
show that we can derive the magnetisation orientation 
of the magnetic contacts and their polarisations. Com- 
parison measurements on a different device at 4.2 K are 
shown at the end of this section. We will also discuss the 
relationship between tunnel barriers' transparencies and 
polarisation. 

To characterise the device, each individual tunnel bar- 
rier is first measured: with reference to Fig.^ by sending 
a current between Col and Co2 and detecting the volt- 
age between Col and Co3, we measure tunnel barrier 1. 
Usually Col and Co3 have the same resistance (within 
20%) as do Co2 and Co4, because they have nominally 
the same area. 

The I-V characteristic of a 1.8 kfl tunnel barrier mea- 
sured at liquid helium temperature is shown Fig.^] Pos- 
itive voltage means Co at higher potential than Al. All 
the tunnel barriers (TB) we measure at 4.2 K (with re- 
sistances down to 7 kn) show a peak at zero bias and 
are asymmetric in the applied bias. Variation of the tun- 
nel barrier differential resistance of 10% in the bias range 
used is visible. 

The measurements are organised in the following way: 
for each measuring configuration, we show one spin valve 
measurement and three precession measurements for the 
different magnetic configurations. 

Figure EJa) shows measurements at RT for the side 
configuration. Starting with the magnetic field at 80 niT, 
with all the magnetic contacts pointing parallel to each 
other, we sweep the field to negative values. At —25 mT, 
the two larger electrodes flip, the magnetic configuration 
is antiparallel and the detected signal increases above the 
background level by 45 mf2. Increasing the field further, 
at —32 mT, one of the smaller electrodes reverses, and 
the signal dips 10 mfi below the background level. At 
—38 mT, the second smallest electrode also flips and the 
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FIG. 5: (a) Spin valve measurement of i? = [V^ — V~)/I as a 
function of the in-plane magnetic field at RT. (b) Memory ef- 
fect starting from the antiparallel configuration (upper trace 
beginning at 246 mf2) and from the anomalous configuration 
(lower trace from 186 mfi). (c) Precession measurements as 
a function of the magnetic field applied perpendicular to the 
substrate for the antiparallel, parallel and anomalous config- 
urations. The fits for the antiparallel and anomalous cases 
follow closely the experimental data. The right inset shows 
the direction of the electrodes' magnetisations. The left inset 
shows precession measurements with the external field applied 
perpendicular to the spin accumulation in the x direction. At 
fields as low as 30 mT, the contacts' domains begin to turn, 
and the signal becomes irregular. 



signal reaches the background value. The reverse trace 
shows very similar behaviour, the notable difference be- 
ing a larger peak, about 48 mf2, and a smaller dip, 6 inQ. 
Repeated sweeps give similar results. 

Figure |3Jb) shows the memory effect, that reflects the 
hysteretic behaviour of the electrodes. Starting with the 
system in the antiparallel configuration at -1-30 mT, we 
sweep the field toward negative magnetic fields. The elec- 
trodes stay in the antiparallel magnetic configuration un- 
til we reach —25 mT, at which point the largest electrodes 
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switch parallel to the smallest ones. In the reverse sweep, 
at +20 mT the largest electrodes flip again, returning the 
initial configuration. The second trace is the memory ef- 
fect in the anomalous configuration. Suppose Col, Co2 
and Co3 are parallel to the y direction and Co4 is oppo- 
site. Starting at -1-33 mT (and 186 mil), we sweep the 
field towards negative fields. At —25 mT, the largest elec- 
trodes Col and Co3 flip, now pointing in the — y direction 
and parallel to Co4: this is still an anomalous configura- 
tion. Upon reaching 40 mT, the smallest electrode Co3 
flips parallel to the other three and the detected signal 
reaches the background level. 

In a precession experiment, we apply an external mag- 
netic field perpendicular to the sample (in the z direc- 
tion), with the in-plane field switched off. Figure Etc) 
shows precession measurements for the three magnetic 
cases, parallel, antiparallel and anomalous configura- 
tions. In both antiparallel and anomalous there is a no- 
ticeable dependence of the signal on the magnetic field, 
whereas in the parallel case little modulation is seen. 
From the smoothness of the curve, we conclude that the 
contacts' domains do not flip irreversibly in the direction 
of the external field, up to fields of 280 mT. At higher 
fields, however, the magnetisation of the end domains of 
the strip is unstable and tends to fiip to a different con- 
figuration irreversibly, resulting in sudden jumps of the 
signal. One could also think of doing precession measure- 
ments with an in-plane magnetic field, perpendicular to 
the leads, in the x direction. The left inset in Fig. [SJc) 
shows the result of such a measurement, a sweep from 
+40 mT to -60 mT and back to -1-60 mT: at fields as low 
as |30| mT, the electrodes' domains begin to rotate and 
the signal deviates from a smooth curve. 

The signal is fit with eq. H16|) , written in a way suitable 
for interpolation, and similar to eq. (2) of Ref. 
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The background term Rback accounts not only for the 
Ohmic resistance, but also for the magnetisation that 
is injected parallel to the applied field and that does not 
precess. In fact, from the geometry of the device, the end 
domain of the contacts are pointing slightly upwards. 

Note that only in the case Wst,_L = 0, is the angle 
between injector and detector (it is actually the projec- 
tion of the angle on the plane perpendicular to z), and 
Trei = Taf. Here, we allow for a stray field through the 
sample only in z direction: lOz = .g/i_B|B|/?i -f Wst^z and 
we assume f~; tJ^. 

For the antiparallel case of Fig. ^p) we find r^/ = 
62 ± 2 ps = (-0.06 ± 0.01)7r, ujst = -14 ± 4 mT 
and a background Rback = 192 ± 1 mfi. Using the dif- 
fusion constant for the aluminium Z? = 5 x 10~^m^/s 
found from resistivity measurements, the diffusion time 
L^/D « 30 ps, shorter than r^/. The spin 
y^Z^rTf = 550 nm. 






FIG. 6: (a) Orientation of the electrodes' magnetisation 
in the parallel configuration. The quantity Pva.± for the 
widest electrodes is canted by an angle (j)' with the direc- 
tion of the narrow electrodes' magnetisation and is shorter 
by K than the narrowest electrodes. The black arrows rep- 
resent the injector sx = (Pi mi — P2ni2)± and the detector 
dx = (Pjma — P4m4)x, <^ being the angle between the two 
vectors sx and dx. The same schematics for the antiparallel 
(b) and anomalous (c) configurations. 



We now give a first estimate of the polarisation P: 
assuming that only the widest electrodes' magnetisations 
are rotated by 0, |s| = |d| = 2Pcos((/)/2), we find P — 
7%. 

Parallel and anomalous configurations differ by about 
6 — 10 mfi in the spin valve measurement, Fig|SJa), and 
12 ToD, in the precession trace, whereas the spin current 
accounts for only P^Rohm — 1 mil. The precession data 
for the anomalous configuration indicate that accumula- 
tion also occurs2i, and accounts for most of the signal. 

We fit the signal in the anomalous configuration with 
eq. (|16|l . but now fixing r^f to the value found in the 
(side) antiparallel case. We find = (0.11±0.04)7r, utst = 
± 8 mT and Rback = 197 ± 1 mfl. 

From the precession measurements, we work out the 
magnetic configuration of each electrode. We note first 
that for the function 
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IS Tdiff 

diffusion length is As 



max((7) — min(g) — 1 holds, for every value of 0: the 
amplitude of the spin signal in a precession measurement 
is proportional to |sx||dj^|, independent of the angle be- 
tween injector and detector. 

We show now that the precession measurements for 
the three magnetic configurations are consistent if one 
assumes that the narrow electrodes Co2 and Co4 point 
in the y direction, that Col and Co3 are tilted inwards 
by an angle (j)' = (0.08 ± 0.03)7r and that their compo- 
nent on the x — y plane, |Pmx| (i.e. the component that 
precesses) is smaller than the narrow ones by a factor 
K = 0.7 ± 0.1, see Fig. El In fact, |si||dj_| « (l-f k)^ 
in the antiparallel case (cos((/)') ~ 1) and w (1 — k^) for 
the anomalous configuration. Their ratio is {1 — k^)/{1 + 
nf (0.3/1.7) = 18%. We now show that this value 
is close to the experimental result. From the measure- 
ments of Fig. EJc) , we find the maximum modulation of 
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FIG. 7: Spin valve measurement (a) and precession measure- 
ments for the three different magnetic configurations (b) in 
the diagonal case. The fit for the antiparallel and anomalous 
configurations are superimposed onto the experimental data. 



FIG. 8: (a) Spin valve measurement and (b) precession mea- 
surements for the three different magnetic configurations in 
the opposite measuring configuration. The anomalous con- 
figuration for the precession measurement is set by applying 
—40 mT in the y direction. 



the precession signal in the anomalous and antiparallel 
configurations, respectively 9 and 55 mQ; their ratio be- 
ing 16%. 

This is also compatible with the small magnetic field 
dependence of the parallel configuration signal. In fact, 
the ratio between magnetic signals in the parallel and 
antiparallel cases (3 mft/db mil w 6%) is close to the 
expected ratio, (1 - + n)^ = (0.30/1.70)2 = 3% 

(cos</)' ~ 1). With these corrections, the efficiency of the 
narrowest electrodes becomes P = 8%. 

Taking into account the efficiency of the detector P, 
the spin accumulation is Rs/P = 850 mil, larger than 
the Ohmic background resistance. With a typical driving 
current of 10 /iA, we find the imbalance between up and 
down spins in the island An = RJiyoeV/P = 10^. For 
comparison, the total number of free electrons is 10^. 

The spin signal in the side configuration is compared to 
that in the diagonal configuration. Fig. [JJa). The signal 
of the spin valve measurement, 54 mfi for the left peak 
and 48 laQ for the right one, is comparable to that in 
the side configuration, supporting the assumption of the 
spin accumulation in the island being uniform. 

Using this measuring configuration, we also perform 
precession measurements in the parallel, antiparallel and 
anomalous magnetic configurations. As before, we find 
the relevant parameters, P = (8.0±0.5)% and Tgf = 65± 
4 ps, consistent with those found in the side configuration 
{Rback = -8±lmn, ujst = -13±4 mT and = (0.10 ± 
0.02)7r). In the anomalous diagonal configuration, we 



fix the relaxation time found in the diagonal antiparallel 
case and we find (/) = 0.08 ± 0.05, uJst = -8 ± 8 mT, 
Rback = -7 ± 1 mil. 

Both LOst and (p show variations between successive 
measurements, spanning from |7| to |15| mT and from 
|0.04|7r to |0.07|7r respectively, in the antiparallel side 
configuration. On the other hand, both the spin relax- 
ation time and the polarisation showed constant values 
throughout the time of the measurements. 

In the opposite configuration, the signal in the parallel 
and antiparallel configurations differ by 5 — 8 mi7. The 
precession measurements show dependence on the exter- 
nal B field of less than 3 mil and we conclude that the 
magnetisation is injected parallel to the external field. 

The signal in the anomalous opposite configuration 
shows a magnetic field dependence. The most notable 
feature is that the signal is odd in the magnetic field, 
implying that s and d are almost perpendicular to each 
other. We fit the signal by fixing the spin relaxation 
time Tsf = 62 ps as found from the side configuration 
and find = (-0.35 ± 0.01)7r, uj^t = -23 ± 7 mT and 
Rback = 200 ± 1 m^!. 

We now compare the predicted signal in the side an- 
tiparallel configuration, (I-I-k)^ and in the anomalous op- 
posite configuration 2 X 2k sin (/)', {1 + k)'^ / Ak sm (j>' = 4.1, 
comparable to the amplitudes ratio of the measured sig- 
nal 55 mfl/lO mn = 5.5. 

We also test the prediction of the reciprocity theorem. 
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FIG. 9: Precession measurements in the antiparallel side con- 
figuration to show the reciprocity theorem: curve (1) is mea- 
sured with magnetic field applied in the z direction, (2) is 
measured after interchange of current and voltage probes and 
magnetic field applied in the — z direction, curve (3) after 
reversal of all magnetisations, interchange of voltage and cur- 
rent probes and magnetic field in the — z direction. 



which states that a four-probe measurement is invariant 
upon exchange of the voltage and current probes and 
magnetic field reversal^-'^. In the case of magnetic elec- 
trodes, one has also to reverse their magnetisations. Fig- 
ure 13 shows measurements of the spin accumulation for 
different electric and magnetic configurations. We pro- 
ceed as follows: with the external field applied in the 
positive y direction, we set the contacts in the antiparal- 
lel configuration. We then measure the precession signal 
(curve 1). Next, we exchange the current and voltage 
probes and we repeat the measurement (curve 2), this 
time applying the external magnetic field in the — z di- 
rection. Next, with the leads interchanged, we apply the 
the field in the — y direction and set the device's mag- 
netic configuration to antiparallel. We then measure the 
precession signal, again with the external field in the — z 
direction (curve 3). We see that curves (1) and (3) are 
identical. Curve (2) is shifted in magnetic field and its 
maximum is 2 mil higher than the other two curves. In 
fact, in the presence of a ^ 0, the last term in eq. Hl()|) 
is not invariant if we simply exchange s <-> d without 
flipping their directions, s ^ — s and d — d. A small 
in-plane stray field of « 20 mT is enough to account for 
a difference of 2 mi7. 

We also note that the spin valve traces for all magnetic 
configurations at zero magnetic field are offset from one 
another by about 1-2 mfl. We believe that again the 
last term of eq. H16() is responsible, as reversing the elec- 

— nij, also causes uj,,t to 



trodes' magnetisations nij - 
reverses. We can exclude a Hall effect generated by the 
leads' magnetic field. In fact, taking the Hall resistance 
Rh = -3.5 • 10"" n ■ m/T for Al in the low field limiliS^, 
the Hall contribution would be RnB/d w 0.1 mil, us- 
ing the thickness of the island d = 30 nm and a field of 
100 mT, which is too small to explain the difference. 
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FIG. 10; Spin valve measurement at 4.2 K in the (a) side, (b) 
diagonal and (c) opposite configurations. The reversal of the 
widest electrodes is step-like. In the diagonal configuration, 
the left and right peak have different height, probably due to 
incomplete switching of one of the contacts. 



As a comparison. Fig. H10|) shows measurements at 
4.2 K on a different device with tunnel barriers of 2 kfi 
(Al deposition rate 0.2 nm/sec). Side and diagonal con- 
figurations show a similar spin signal, around 230 mfi. 
In the diagonal case, the left peak does not reach full 
height, probably due to the incomplete reversal of one 
of the wide electrodes. The switching of the magnetisa- 
tion occurs with discrete changes, resulting in a step-like 
spin signal, as opposed to the switching at RT, which oc- 
curs abruptly in most cases. The opposite configuration 
shows little spin signal, as expected. 

We have also measured two devices with tunnel bar- 
riers in the range 2-4 kO, and found a spin signal of 
80 mil. For more transparent interfaces, 0.8-1.6 kO (two 
devices), the spin signal is 55 mfi. For the last four de- 
vices, the spin relaxation time is Tsf = 60 ± 4 ps. Two 
devices (Al deposition rate 0.2 nm/sec) were measured 
both at 4.2 K and at RT: one device, with tunnel barri- 
ers of 5-11 kfi, gave a spin signal of 90 mfl at RT and 
250 mfi at 4.2 K, the other with tunnel barriers 15-35 kil, 
150 mfl at RT and 300 mfl at 4.2 K. 



12 



We note that the polarisation of the interface decreases 
from 10.5 % for the highest resistance interfaces to 8 % 
(if we assume that the relaxation time is the same for all 
devices). Although these values are unequivocally lower 
than spin polarisation measurements with superconduct- 
ing aluminium:3&, we see that the aluminium oxide inter- 
faces can be made transparent enough without losing the 
polarisation completely. 

Jedema et ali has performed for AI/AI2O3/C0 1-D 
structures, spin valve and precession measurements. E- 
beam evaporation was used (the same evaporating ma- 
chine) to deposit the metals. They also reported spin flip 
times of 50 ps but a polarisation of the tunnel barriers 
at RT of 11 %, slightly larger that what we found, based 
on the fit to the experimental traces with a time-of-flight 
1-D model that takes the diffusion constant as an in- 
dependent parameter. They also found, from precession 
measurements at 4.2 K, that the spin polarisation of the 
tunnel barriers increases to 13 % and the spin relaxation 
time doubles. 

As we were not able to perform precession measure- 
ments at 4.2 K (due to technical difhculties), we could 
not determine separately the value of the spin relaxation 
time and the spin polarisation of the injector and detec- 
tor. For this reason, we cannot conclude which is the 
main mechanism of relaxation, whether phonon or im- 
purity induced, and to compare with the theoretical cal- 
culations of Fabian and Das Sarma™. Nevertheless, our 
spin valve risults at 4.2 K are consistents with those of 
Jedema et al. 



VII. CONCLUSION 

Spin accumulation is analysed for zero-dimensional 
systems, in which the electron spin diffusion time Tdiff 
is shorter than the spin relaxation time Trei and the spin 
accumulation can be considered uniform. In the system 
under study, spins are injected into a small island of nor- 
mal metal through ferromagnetic contacts, and the re- 
sulting magnetisation is electrically detected by means of 
other FM contacts. We have theoretically modelled the 
island using the finite element theory of Brataas et aLU.: 
we have shown that the presence of the leads affect the 
spin accumulation by making available extra channels of 



spin relaxation. In particular, the mixing term G'^^ is se- 
lectively relaxing spins with orientation perpendicular to 
the electrode magnetisation. The expression we derived 
for the spin accumulation in the island is valid in the case 
of negligible spin accumulation in the FM contacts. 

Experimentally, we have fabricated a small island of 
aluminium with all dimensions (400 nmx400 nmx30 nm) 
smaller than the spin relaxation length (A^/ = 550 nm at 
RT). Transparent tunnel barriers between the island and 
the FM electrodes provide a spin dependent resistance 
that is much higher than all the other (spin indepen- 
dent) resistances in the system, so as to overcome the 
conductivity mismatch. Because of the lateral dimen- 
sions of the island compared to the spin flip length, only 
pure spin accumulation occurs in our device: the spin 
signal can therefore be described in terms of the relative 
orientations of the magnetic electrodes. Spin valve and 
spin precession measurements were presented for differ- 
ent electrical configurations. The peculiarity in the ex- 
periments is that the Ohmic drop across the island is 
smaller than the spin signal. 

In spin valve and precession measurements, we extract 
the polarisation of the tunnel barriers and the spin flip 
time. The AI2O3 tunnel barriers, with resistances of 20 — 
100 ilfim.^ still present a certain degree of polarisation, 
P — 8%. The spin relaxation time at room temperature 
was found to be Tgf = 60 ps. 

The presence of tunnel barriers confine the electrons 
inside the island and they tunnel out of the system long 
after having lost their spin information: in fact the escape 
time is three orders of magnitude larger than the spin 
flip time. With our lateral devices, it is not possible to 
directly measure the spin mixing conductance. 
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